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^(The  effects  of  a  tenuous  exterior  plasma  corona  with  anomalous 
resistivity  on  the  compression  and  heating  of  a  hollow,  collisional 
aluminum  z-pinch  plasma  are  predicted  by  a  one-dimensional  code.  As  the 
interior  ("core")  plasma  is  imploded  by  its  axial  current,  the  energy 
exchange  between  core  and  corona  determines  the  current  partition.  Under 
conditions  of  rapid  core  heating  and  compression,  the  increase  in  coronal 
current  provides  a  trade-off  between  radial  acceleration  and  compression, 
which  reduces  the  implosion  forces  and  softens  the  pinch.  Combined  with 
an  heuristic  account  of  energy  and  momentum  transport  in  the  strongly 

c 

coupled  core  plasma  (typically  0.1)  and  an  approximate  radiative 

loss  calculation  including  A1  line,  recombination  and  Bremsstrahlung 
emission,  the  current  model  can  provide  a  reasonably  accurate  description 
of  imploding  annular  plasma  loads  that  remain  azimuthal ly  symnetric.  The 
implications  for  optimization  of  generator  load  coupling  are  examined.  { 
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(e)  Peak  Compression 


•  T,  ~  0.5  keV, n-j  ~  1020  cm'3 


Reversal  of  Inward  Momentum 


•  Peak  of  l2,  T2 


(f)  Late  Time  Behavior 


Simple 


— -- rrrrrr:.  ^  ;.  •  Hot  Spot 

ggl*  (Sausage  Instability) 

•  Ribbon 
(Kink  Instability) 


•  Requires  Most  Detailed  Treatment  of 
Plasma  Dynamics  and  Radiation 


Dynamics  in  the  Si 


Model 


Energy  Out 


Radiation  Losses 


Radial  KE, 
Momentum 
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Generator 

V„(t) 


Effective  Circuit  Relation 


The  Circuit  Element  Replacing  Vp)  is  Like 


With  the  Current  Source  l2(t)  Controlling  the  Inductance  and  Resistance  Values. 


II.  COROflA  EVOLUTION: 


Assume  a  marginally  stable,  drift-speed  limiting  system  in 

AN  AVERAGED  STATE  OF  PRESSURE  BALANCE  - 


v(vn[V  <TJ)  =  c''J2(r)XBdCr) 

|4sfr)|-%(0*e/O 


FOR  THE  ISOTHERMAL  SYSTEM,  T 2  AND  JZ(A)  DETERMINE 

{  ne  (r) ,  I  =  !,  +  !,  } 


UNIQUELY. 


WHERE 


03  Y"  =  [l  -0.  +  (30o/2)'] 

ne(r)=  n  /tY-*ri±4 _ 

aU/  \(r/*)1  +  S 


~  3P./2 


Ii/ewCj  (a'-rf) 

3  3  (n/»)7fi 

A>  2  zc^Tj  /ec^I  . 

IF  I(t),  g(t)  ARE  GIVEN,  ONE  NEEDS  ONLY 

jf(N2Z2  Tj  =  1 0  */s)Ez  +  Q|Z  -  Q„a  -  Qen<l 
3?  (Nt )  =  F.z  *  Z,"‘ 


(EFFECTIVE  ION  EXCHANGE  RATE) 
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Qn  =  5"™vt  T« 


F  =  Wav, 

i  * 


WITH  Fj2  evaluated  by  means  of  a  simple  Krook  model  in 
THE  TWO  REGIONS-. 

n*Tev 


Z  m. 


fleTe 


2"*m« 


[  (-drMnt) 

T«w  Cdr -fcn 

l  1  *  (n.Tee)1 

1  <-  / 

(  T«(-3ri»tfie) 

\  1  ♦ 

l  +  (n„  T'J1  1 

i 


IN  ORDER  TO  ADVANCE  THE  COMPLETE  CORONA  IMPLICITLY  FROM  ONE 
PRESSURE  BALANCE  TO  THE  NEXT* 


i 

i 
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III.  CORE  PLASMA  EVOLUTION 


A.  Strongly-coupled  plasma 

SpITZER's  PARAMETER  Mi  A  ,*j  =  23  - Xn.\[znL  Z3(Tc }/Tj3/2'] 
(I;,  Te  in  €.v  ,  r>.  in  c.m’5  j 

•  Initially  15  eV,  z  ~3,  3  x  1018^ni^1.2  x  109 

MiA n  ~  2.H  J.Hf 

•  Typical  peak  compression:  750  eV,  z  «  10, 

(l-^ij)x  1019  CM-3  : 

Mi  A  n  ^  3.8  —  3.1 

20 

§  Later  collapse  phase:  30 eV,  zs:  4,  (Mjx  10  : 
sU  A  n  x  o.t>  -*  -  o.i2 

Classical  transport  applies  to  plasma  with  -lu-Au^lO* 

Alternative:  dissipative  processes  dominate  fluid 
evolution: 

(i)  RAPID  THERMAL  CONDUCTION  - *>  V“JT  ^  0 

( 1 1 )  RAPID  THERMAL  EQUILIBRATION - ►  Te  T; 

(ill)  RAPID  VISCOUS  MOMENTUM  TRANSFER - *•  v(v»v);=sO 

(IV)  RAPID  RADIATIVE  DIFFUSION - »->  v(n;#/ rvJ^O 

Result:  simplified  radiation/hydro, 

USEFUL  MODEL  WHERE  CLASSICAL  THEORY 
INAPPLICABLE 


B.  Initial  Conditions: 


•  Choose  compatible  with  (i-iv). 

•  Current  in  outer  part  of  core:  (*i,a.) 

•  Vn;/Di  oc  Jzx  Be/ne(r)  -  otv(r) 

•  V(r)  =  Gj  r  *  C-,_J  r"1 

(These  profiles  remain  self-similar  in  slab-geometry 
expansion ,  compression  and  translation,  provided 
Ez(r,t)  holds  Jz/cne-v,  constant.) 

Parameter  cc  controls  DVr/Dt  at  *t0 . 

Evolution  without  external  stress:  Dvr/Dt=  to  give 

viscous  stress  relaxed  to 
zero. 

(Density  maximum  stationary  in  this  case,  while  pro¬ 
file  EXPANDS  SUBJECT  TO  BOUNDARY  CONDITIONS  AT 
a  2,  b  -  WHERE  Vp  MOVES  PLASMA.) 

C.  Response  to  External  Forces 

•  Transmission  of  stress  from  current-carrying  layer 

to  interior  fluid:  f>£r  +  c-  JzX  Be  +  . 

F  arise  in  maintaining  v*v  homogeneous 

Tj  Vise 


•  Classical  plasma  can  snowplow;  strongly-coupled  viscous 

PLASMA  ERASES  INHOMOGENEOUS  COMPRESSIONS  AND  ABSORBS 


MOMENTUM  OVER  ENTIRE  VOLUME. 


•  Simplest  model:  outer-surface  stresses  are  transmitted 

UNATTENUATED  INTO  PLASMA.  VlSCOUS  REACTION  FORCES 
'INSTANTANEOUSLY'  REDISTRIBUTE  MOMENTUM  TO  KEEP  v(V*v)  =  Q 


•  Simple  representation  of  model:  focus  on  density 

MAX  (rm)  AND  INNER  RADIUS  G>) ,  WHERE  EFFECTIVE  STRESSES 
ARE  KNOWN: 


^  la.  ) 


P  =  -  X  ¥IL 

b  m  O  b 

•  Homogeneity  of  v»v  then  requires 


AND  GIVES  T>Vr/Pt  =:  <*i(v.-rm)  +  o(0 


D.  Assembly  on  Axis: 

Strong  vn  ,  vv 

Interpenetrating  strongly-coupled  fluids 

Microturbulence 

- > 

Fluid  mixing 

Dissolution  of  vn  i-  vv 

Thermalization  of  flow  kinetic  energy 


m 


Use  simple  "collection"  model: 

(1)  Inside  a  "collection  radius"  rc  ,  thermalize  any 

FLOW  ENERGY  IN  EXCESS  OF  THAT  IN  STAGNANT  FLOW 
(V,  =£«  ) 

(2)  When  particles  collected  within  rt  are  numerous 

ENOUGH  TO  ERASE  OR  REVERSE  Vn  AT  Tc  ,  MOVE  rc 
OUTWARD  ONE  GRID  CELL. 

(3)  Continue  doing  (2)  until  rm  meets  rc 

(4)  Continue  compression  of  uniformly-filled  center 

( r<  rc ) 


E.  Motion  of  the  Collected  Plasma 

•  Simple  representation  of  stress  transmission 
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T  *  B  I 

~ (explicit  pressure  difference) 
cmne(rfn) 

( v*v  homogeneity  ) 


§  Self-similarity  of  v(r)  maintained  by 
Dv(r) 


Pi 


-  UV(r) 


FROM  VISCOUS  RESPONSES  RATHER  THAN  FROM  A  PRESCRIBED 


PRESSURE  GRADIENT. 


F.  Core  Energy  Transport  (spatially  integrated) 


Z(Te)  ASSUMED 


Jt  (ZJTl)  =  37^1  *  Qvi«  +  Q 


1  Amy 

STr(a*~  bx)f 


WITH  Q12  THE  CORE/CORONA  INTERFACE  HEAT  EXCHANGE, 

g  • 

QVISC  CALCULATED  FOR  THE  VELOCITY  FIELD  VR(*>, *,b), 

PrAd  a  COMBINATION  OF  LINE,  RECOMBINATION  AND  BREMSSTRAH  - 
LUNG. 

G.  External  Circuit 

E z  =  VpLtewX* IN  THE  C0R0NA  RE6I0N 

'Ij.  J(*)-  E(a)+ c  ON  core/corona  interface. 


From  Faraday's  law 
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Given  I(t),  V^Ct),  a(f),  ^(t),  and  the  corona  parameters 

THE  EXPLICIT  DERIVATIVE  I (t)  IS  CALCULATED. 

H.  Tracking  the  density  profile 

Let  Vft)2  iTi.[^(t)  -  b*(t)]  ,  AND  IF 

v(v-v)  =  o  i-  r.  =  rft)  +  f  Vtt  v(r(i,)}  t, ) 

THEN 

rt 

•  n(r(t),  t)  -  n(r0,t0)exp-J  dt^v.vjfr^tt),  tj] 

or  nfr.t)  _  Ylffo) 
no(r.,i *)  ~  Vtt)  * 

•  ANY  n(r«,  )  CAN  BE  EVOLVED  BY  FOLLOWING  r(t,r0). 


AND  INCLUDING  THE  BALANCE  OF  THE  CIRCUIT 

Vt,s  iz9  IR'  - 

WITH 

-Mn-faLiY1)  ~  ^r)  . 

fi,  o t,  6  CORONA  PARAMETERS. 


IV.  RADIATIVE  LOSSES 


For  diagnostics.,  divide  into 

•  Low  ENERGY  LINES  (  h^<lMV) 

•  High  energy  lines  (hz/>lkeV  ) 

•  Continuum  (recomb,  and  Bremss.) 

Use  curve  fit  to  cre  model  of  Duston  and  Davis*1^  for 
A!.  (Thin  cylinder ,  .05  cm  radius  elemental  radiator.) 

Use  Apruzese  probability-of-escape  to  further  attenu¬ 
ate  EACH  SPECTRAL  CATEGORY. 

Integrate  radiative  loss  over  the  density  profile. 


V.  TYPICAL  IMPLOSION  TRAJECTORIES 


•  Simple  Collapse  (SC)— implosion  heating  too  weak  to 
reverse  radial  motion.  Density  peak  after  temperature 

PEAK. 

iJEL 


•  Bounce  and  Collapse  (BO— implosion  heating  halts 

RADIAL  MOTION.  PLASMA  BOUNCES.,  COOLS,  COLLAPSES  AGAIN. 

Temperature  and  density  peaks  closer  together 


0R_ 


•  Pause  and  Collapse  (PC)— on  the  dividing  line  between 
SC  and  BC.  Plateau  in  density,  then  further  compres¬ 
sion. 

•  Bounce  and  Dissipate— Not  seen  for  the  driving  voltage 

WAVEFORMS  USED.  RADIATIVE  LOSSES  PREVENT  THIS 
TRAJECTORY. 

§  The  typical  initial  condition  is  shown-  •'  -••• 
rm»~0.?6 
9  1-0 
71  (o)  ~  Iffe-V 
7Sx  10“  A l  /TOMS 

I* 

•  SC  and  BC  Implosions:  Figure  shows  Lagrangian  zones 

[b(t)  <  rj  ft) <  rm( t)  <  <  aft)  ] 

AVERAGE  DENSITY 

njft) 

ION  TEMPERATURE 

T;W 

As  Functions  of 

-fc  IN  ARBITRARY  UNITS,  VaTH  EARLY  DENSITY  VALUES  SHOWN 
FOR  EACH  ZONE, 

•  A  SEQUENCE  OF  BC  DENSITY  PROFILES,  (N)  DENOTES  THE 


NUMBER  OF  TIME  STEPS. 


Pr-  (tf 9/cm3  mc)  M  a  function  of  T#  (ktV) 

Ifrom  a  0.06  cm  radiut  cylindrical  Al  plawna  with  n.  ■  1.0M01* 


f* 

T 

S 


Tj[eV] 


Radial  Profiles  of  a  Bounce  and  Collapse 
Implosion  [xxx]  Time  Sequence 


[1601 


VI.  GENERATOR  AND  LOAD  COUPLING 


As  A  GENERAL  MEASURE  OF  IMPLOSION  PERFORMANCE;  DEFINE 
T  2  SUP  7;(t) 

OVER  THE  TIME  SCALE  tMAX  CHARACTERISTIC  OF  A  SINGLE  COMPRES¬ 
SION  AND  SUBSEQUENT  EXPANSION  OF  THE  LOAD. 

IN  THE  FIGURE  BELOW  THIS  PARAMETER  IS  SHOWN  AS  A  FUNCTION 

OF  .* 

R°  -  INITIAL  DENSITY  PROFILE  PEAK  (LOAD  RADIUS) 

Vgen  "  PEAI<  GENERATOR  VOLTAGE 

5°  -  INITIAL  CURRENT  PENETRATION  PARAMETER 

Mq  "  TOTAL  MASS  OF  THE  LOAD 

AND  THE  COMMON  CASE  AMONG  THESE  FIGURES  IS 

ss  0.96  (cm) 

=  VQ  (ARBITRARY  units) 

=  9.0 

=  12,  1.5  mil  AX  wires  (7.51  x  1018  atoms  of  A£) 

The  absolute  location  of  these  curves  is  sensitive  function 

OF  THE  CURRENT  PENETRATION  PARAMETER  3°,  WHICH  EFFECTIVELY 
SETS  A  "MEAN  SKIN  DEPTH"  FOR  THE  ELECTROMAGNETIC  FIELD  PENETRA¬ 
TION  OVER  THE  IMPLOSION  HISTORY. 


T,  teVD 


m  t arb . ] 

m=l  is  1.SE19  fil  atoms 


VII.  CONCLUSIONS 


The  implosion  dynamics  are  dominated  by  the  amount  of  energy 

DELIVERED  BY_ THE  GENERATOR  TO  THE  RADIAL  FLOW  OF  THE  FLUID  LOAD. 

•  For  low  values  of  this  energy,,  the  heating  on  assembly  is 

TOO  WEAK  TO  REVERSE  THE  IMPLOSION.,  AND  A  REFRIGERATION  (SC) 
OCCURS,  WITH  THE  PEAK  COMPRESSION  AFTER  THE  PEAK  TEMPERATURE. 

•  AS  DELIVERED  ENERGY  IS  INCREASED,  THE  COMPRESSION  PEAK  AND 
ASSEMBLY-HEATING  TEMPERATURE  PEAK  BECOME  CLOSER  IN  TIME. 

Assembly  heating  and  compression  heating  work  together  to 

PRODUCE  HIGHER  PEAK  TEMPERATURE. 

•  In  CONTRAST,  WHEN  FURTHER  FLOW  ENERGY  IS  GIVEN  TO  THE  LOAD, 

THE  ASSEMBLY  HEATING  TENDS  TO  HALT  THE  COMPRESSION  VERY 
RAPIDLY  (AND  AT  A  LARGER  RADIUS)  AND  THE  COMPRESSION  REVERSES 

too  soon.  This  results  in  a  lower  peak  temperature. 

All  these  effects  are  superimposed  on  a  very  sensitive 

COUPLING  TO  THE  CORONAL  REGION,  AS  MEASURED  IN  PART  BY  THE 
RELATIVE  SKIN  DEPTH  PARAMETER  3 .  As  THE  CURRENT  IS  FORCED  TO 
A  THIN  LAYER  (g^»)  THE  SURFACE  STRESSES  ON  THE  CORE  PLASMA 
BECOME  STRONGER  UNTIL  VERY  LARGE  3  VALUES  DEMAND  THAT  fi-*0  AND 
FORCE  MOST  OF  THE  CURRENT  TO  THE  CORONA.  As  3-0  THE  CORONA 
BECOMES  RELATIVELY  UNIMPORTANT  AND  THE  SURFACE  STRESSES  WEAKEN. 

An  accurate  assessment  of  the  importance  of  turbulent  current 

TRANSPORT  IN  THE  LOW  DENSITY  REGIONS  OF  THE  ANNULAR  LOAD  MUST 
THUS  DEPEND  UPON  A  PROPER  ELECTROMAGNETIC  TREATMENT  OF  THE 
E,  E. AND  FIELD(^  PROFILES  WITHIN  THE  LOAD,  THIS  IS  A  SUBJECT 
OF  FURTHER  WORK  AND  THE  OUTCOME  WILL  PROVIDE  A  FIRST-PRINCIPLES 
GUIDE  TO  THE  APPROPRIATE  3  VALUES  FOR  THIS  SIMPLE  MODEL. 


REFERENCES 


1.  D.  Duston  and  J.  Davis.,  Phys.  Rev.  A,  2L  Hay  1930. 

2.  J.  Apruzese  and  J.  Davis,  "Direct  Solution  of  the  Equa 
tion  of  Transfer  Using  Frequency-and-Angle -Averaged 
Photon  Escape  Probabilities."  NPL  Hemo  Peport  Number 
4017,  21  June  1979. 

3.  R.  E.  Terry  and  J.  U.  Guillory,  JAYCOR  Report  Number 
TPD200-80-012,  October  1980. 


